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' 1 Introduction 

' Heavy quarkonium systems have been a good laboratory 
] to test our theoretical ideas since the early days of QCD 

■ (see [1] for an extensive review). Based on the pioneering 
' work of Caswell and Lepage [2] , a systematic approach to 
, study such systems from QCD has been developed over 

■ the years which makes use of effective field theory tech- 
\ niques and is generically known as Non-Relativistic QCD 
- (NRQCD) \^ . NRQCD has been applied to spectroscopy, 
' decay and production of heavy quarkonium. The NRQCD 
, formalism for spectroscopy and decay is very well under- 
' stood (this is also so for electromagnetic threshold pro- 
duction) so that NRQCD results may be considered QCD 

■ results up to a given order in the expansion parameters 
' (usually aslrnq) and l/mg, mg being the heavy quark 
, mass). I will restrict myself to review recent progress on 

■ these issues and refer the reader to the review [¥ for ear- 
\ lier developments. The NRQCD formalism for production 

is more controversial. Important work has been carried 
out recently to which I will devote some time. 



2 Heavy Quarkonium 

Heavy quarkonia are mesons made out of a heavy quark 
and a heavy antiquark (not necessarily of the same flavor) , 
whose masses are larger than Aqqd, the typical hadronic 
scale. These include bottomonia (66), charmonia (cc). Be 
systems (6c and cb) and would-be toponia (it). Baryons 
made out of two or three heavy quarks share some simi- 
larities with these systems (see [5j|6]). 

In the quarkonium rest frame the heavy quarks move 
slowly (u 1, f being the typical heavy quark velocity 
in the center of mass frame), with a typical momentum 
rnqv <C mg and binding energy ~ raqv^ . Hence any study 
of heavy quarkonium faces a multiscale problem with the 
hierarchies mg ^ rnqv ^ mqv^ and mg 3> Aqcn- The 
use of effective field theories is extremely convenient in 
order to exploit these hierarchies. 



3 Effective Field Theories 

Direct QCD calculations in multiscale problems are ex- 
tremely difficult, no matter if one uses analytic approaches 
(i.e. perturbation theory) or numerical ones (i.e. lattice). 
We may try to construct a simpler theory (the effective 
field theory (EFT)), in which less scales are involved in 
the dynamics and which is equivalent to the fundamental 
theory (QCD) in the particular energy region where heavy 
quarkonia states lie. The clues for the construction are: (i) 
identify the relevant degrees of freedom, (ii) enforce the 
QCD symmetries and (iii) exploit the hierarchy of scales. 
Typically one integrates out the higher energy scales so 
that the lagrangian of the EFT can be organized as a 
series of operators over powers of these scales. Each op- 
erator has a so called matching coefficient in front, which 
encodes the remaining information on the higher energy 
scales. The matching coefficients may be calculated by im- 
posing that a selected set of observables coincide when are 
calculated in the fundamental and in the effective theory. 

NRQCD is the (first) effective theory relevant for heavy 
quarkonium. Out of the four components of the relativis- 
tic Dirac fields describing the heavy quarks (antiquarks) 
only the upper (lower) are relevant for energies lower than 
mq (no pair production is allowed anymore). Hence a two 
component Pauli spinor field is used to describe the quark 
(antiquark). The hierarchy of scales exploited in NRQCD 
is mg 3> mqv,mqv^ , AqcD- The remaining hierarchy 
^qv 3> mqv^ , may be exploited using a further effective 
theory called Potential NRQCD (pNRQCD) [7118]. 

3.1 Non-Relativistic QCD 

The part of the NRQCD lagrangian bilinear on the heavy 
quark fields coincides with the one of Heavy Quark Effec- 
tive Theory (HQET) (see [9] for a review), and reads 

I 2mQ 8mg 2mg 
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ip is a Pauli spinor which annihilates heavy quarks, cp, cd 
and Cs are short distance matching coefficients which de- 
pend on uiQ and fi (factorization scale). Analogous terms 
exist for x, a Pauli spinor field which creates antiquarks. 
Unlike HQET, it also contains four fermion operators, 



element is necessary to cancel the factorization scale de- 
^pendence which arises in one loop calculations of Im/i 
^2S'+ip^-j . The matrix elements cannot be calculated 
in perturbation theory of asimq). They can however be 
calculated on the lattice (see for instance [13]) or extracted 
from data. The imaginary parts of the matching coeffi- 
cients of the dimension 6 and 8 operators are known at 
NLO [14 , and those of the electromagnetic decays of di- 
mension 10 at LO p^lTBllTTl. 



+^4^0sCSo) + ^^OseS,), 



Oi(i5o) =i^HxH , OiCSi)^ij^cTxx^<Tij, (1) 
08(1^0) = V-tT-xX^T^V , OsCS,) = iP^T'^axx^T'^rrij. 

The fs are again short distance matching coefficients which 
depend on mg and which have imaginary parts when 
the quark and the antiquark are of the same flavor. This 
is due to the fact that hard gluons (of energy ^ mg), 
which remain accessible through annihilation of the quark 
and antiquark, have been integrated out. The fact that 
NRQCD is equivalent to QCD at any desired order in 
1/toq and asimq) makes the lack of unitarity innocuous. 
In fact, it is turned into an advantage: it facilitates the 
calculation of inclusive decay rates to light particles. 

The lagrangian above can be used as such for spec- 
troscopy, inclusive decays and electromagnetic threshold 
production of heavy quarkonia. Spectroscopy studies have 
been carried out on the lattice and are discussed in [TUj 
(see also jTIj and references therein). 



3.1.2 Problems 

Unlike HQET, in which each term in the lagrangian ^ 
has a definite size (by assigning ^gp/j to any operator 
of dimension n), NRQCD does not enjoy a homogeneous 
counting. This is due to the fact that the scales mgv, 
mqv^ and Aqcd are still entangled. Even though a count- 
ing was put forward in [3], the so called NRQCD velocity 
counting, which has become the standard bookkeeping for 
NRQCD calculations, there is no guarantee that it holds 
for all the states, and hence there is a problem on how 
to make a sensible organization of the calculation for a 
given state. There are two approaches which aim at disen- 
tangling these scales, and hence to facilitate the counting. 
The first one consist of constructing a further effective 
theory by integrating out energy scales larger mgw^. This 
is the pNRQCD approach first proposed in [7], which will 
be described in the next section. The second one consist 
in decomposing the NRQCD fields in several modes so 
that each one has a homogeneous counting. This is the 
vNRQCD approach first proposed in [18] (see also [19]). 
The scale Aqcd has not been discussed so far in this ap- 
proach, and, in fact, a consistent formulation has only be- 
come available recently [20] (see also [2T1I22] ). 



3.1.1 Inclusive decays 



3.2 Potential NRQCD 



Let us just show, as an example, the NRQCD formula 
for inclusive decays of P-wave states to light hadrons at 
leading order 



{XQinJS)\0,{^'^+^Pj)\XQinJS)) 



+linfsC''+'Ss){xQinJS)\OsCSo)\xQ{nJS)) 



As mentioned above, pNRQCD is the effective theory which 
arises from NRQCD after integrating out energy scales 
largen than mgw^, namely than the typical binding en- 
ergy. If Aqcd < wgw^, then ttiqv > Aqcd and the 
matching between NRQCD and pNRQCD can be car- 
ried out in perturbation theory in as(fTigw) . This is the 
so called weak coupling regime. If Aqcd ^ "T-gw^, the 
matching cannot be carried out in perturbation theory in 
cts{mQv) anymore, but one can still exploit the hierarchy 
mQ ^ mQV, Aqcd 3> mQV^ . This is the so called strong 
coupling regime. 



/i and /s are short distance matching coefficients which 
can be calculated in perturbation theory in as(r7ig). Oi 
^S+ipj-^ is a color singlet dimension 8 operator and Os(}-Sq 
is the color octet operator given in ([T]). Earlier QCD fac- 
torization formulas were missing the color octet contribu- 
tion. They are inconsistent because the color octet matrix 



3.2.1 Weak Coupling Regime 

The lagrangian in this regime reads 



pNRQCD 



= / d^rTr<^ST(«5o-/z,(r,p,PR,Si,S2,M))S- 
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+0t [iDa - Kir, p, Pr, Si, 82,^)) o| 

+VA(r,A*)Tr{Otr -gES + S^r -^EO} + 
+ YEhj^T, {otr . .gE O + OtOr • .9E} 

where 5* and O are color singlet and color octet wave func- 
tion fields respectively, hs and ho are color singlet and 
color octet hamiltonians respectively, which may be ob- 
tained by matching to NRQCD in perturbation theory in 
asijnQv) and l/mg at any order of the multipole expan- 
sion (r ~ 1/mQv). The static potential terms in hg and ho 
are known at two loops p5[|24l[^ (the logarithmic con- 
tributions at three loops and (for kg only) at four loops 
are also known 26,,8.,27j ) . The renormalon singularities in 
the static potentials are also understood in some detail 
[28l[29] . The l/mq and 1/?tiq terms in hg are known at 
two and one loop respectively [30j. Va,Vb = 1 + 0{al) 
[27] . This lagrangian has been used to carry our calcula- 
tions at fixed order in as- an almost complete NNNLO 
(assuming Aqcd ^ 'tiq^s) expression for the spectrum 
is available [31]. Most remarkably resummations of loga- 
rithms can also be carried out using renormalization group 
techniques [3l|33l[34] Q Thus the rjb mass and its electro- 
magnetic decay width have been predicted at NNLL and 
NLL respectively [551155] . 



3.2.2 Strong Coupling Regime 

The lagrangian in this regime reads Q 

ipNRQCD = y^'R J d^'r S\ido-hg{r,p,Pn,Si,S2))S, 



were i?„(0) is the wave function at the origin, Bi ^ ^qcd 
is a universal (independent of n) non-perturbative param- 
eter, and cp a (computable) short distance matching co- 
efficient. This further factorization allows to put forward 
new model independent predictions, for instance the ratios 
of hadronic decay widths of P-wave states in bottomonium 
were predicted from charmonium data [43| , or the ratio of 
photon spectra in radiative decays of vector resonances 
l46l. as we will see below. 



3.2.3 Weak or Strong (or else)? 

Since mg , v and Aqcd are not directly observable, given 
heavy quarkonium state it is not clear to which of the 
above regimes, if to any , it must be assigned to. The 
leading non-perturbative (~ Aqcd) corrections to the 
spectrum in the weak coupling regime scale as a large 
power of the principal quantum number ^47i48j, which 
suggests that only the n = 1 states of bottomonium and 
charmonium may belong to this regime. However if one 
ignores this and proceeds with weak coupling calculations 
one finds, for instance, that renormalon based approaches 
at NNLO [491150] and the NNNLO calculation [36] give a 
reasonable description of the bottomonium spectrum up 
to n = 3. It has recently been proposed that precise mea- 
surements of the photon spectra in radiative decays, as 
the ones carried out by CLEO [51] , will clarify the assign- 
ments [46] . It turns out that in the strong coupling regime 
one can work out the following formula, which holds at 
NLO, 



dTn 

dz 

dFr 

dz 



Cl [^Sl] [Z) TUQ 



{En — Er 



/i,(r,p,PR, 81,82) 



P 



J"R2 

Aran 



Vs{v, p, Pr, 81, 82), (Oi(35i))„ _ r {T{n) 



{0,{^S,))r r{T{r) 



^ Imgee i^Si) En - Er 

Im/ee (^^i) mg 



Vs 



(1) 



Vs 



(2) 



mq 



Vs cannot be calculated in perturbation theory of as{mQv) 
anymore, but can indeed be, and most of the terms have 
been, calculated on (quenched) lattice simulations 37,38] 
[39 ] I40 P 41] (see [42] for a large Nc calculation). The fact 
that Aqcd ^ fnQV^ can now be exploited to further fac- 
torize NRQCD matrix elements [43ll44ll45] . for instance 



(T(n)|08(i5o)|T(n)) =Ca 



|fl»(0)| 
27r 



(C7^/2 - C/)4gi 
3m| 



^ The importance of log resummations was first emphasized 
in the vNRQCD approach ^18 . However, the correct results 
for the spectrum and electromagnetic decay widths were first 
obtained in pNRQCD |32II33| and later on reproduced in vN- 
RQCD [S]. 

^ We ignore for simplicity pseudo-Goldstone bosons. 



C{ {z), Ci P^i] (z), Imgee (^^i), Im/^e i^Si) are match- 

ing coefhcients computable in perturbation theory, and 
En — Er the mass difference between the two states. If 
data follow this formula it will indicate that both n and 
r are in the strong coupling regime. For the n — 1,2,3 
of bottomonia, current data disfavor n = 1 in the strong 
coupling regime and is compatible with n = 2, 3 in it. 



3.3 Soft-Collinear Effective Theory 

I For exclusive decays and for certain kinematical end-points 
of semi- inclusive decays, NRQCD must be supplemented 
with coUinear degrees of freedom. This can be done in the 
effective theory framework of Soft-Collinear Effective The- 
ory (SCET) 52,53J. Exclusive radiative decays of heavy 
quarkonium in SCET have been addressed in [54] . where 



^ States close or above the open flavor threshold are expected 
to belong neither the weak nor the strong coupling regimes. 
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Fig. 1. is the number of photons and z the energy of the 
photons normalized to half the T{1S) mass. The black squares 
are CLEO data [51] and the (red) solid curve the theoretical 
result of [63] with theoretical errors (grey band) displayed. The 
(green) rectangle on the right displays the region where the 
theoretical results are not reliable. The plot is taken from "64^ . 



results analogous to those of traditional light cone fac- 
torization formulas have been obtained [55] . Concerning 
semi-inclusive decays, let me focus on the photon spec- 
trum of T{1S) — > 7X when the photon is very energetic. 
Describing experimental data in this region has been a 
challenge since the early days of QCD, and the best fits 
were obtained using a model with a finite gluon mass of the 
order of the GeV. [SF. The naive inclusion of color octet 
contributions in the NRQCD framework seemed to further 
deepen the discrepancy [SJ- In recent years a remarkable 
improvement in the understanding of this end-point region 
has been achieved by combining SCET and NRQCD. The 
following factorization formula has been proved [58(159] 



— — = ^ iJ(Af, w, /x) / dkS{k,^i)lmJ^{k + M{l~ z),^i) 



where H is a hard matching coefficient computable in per- 
turbation theory, is a jet function (also computable 
in perturbation theory if z is not too close to 1) and S 
is a (soft) shape function which depends on the bound 
state dynamics. M is the heavy quarkonium mass and ji 
a factorization scale. Large (Sudakov) logs have been re- 
summed using renormalization group equations both in 
the color octet [6^ and color singlet [BT contributions. 
Furthermore, using pNRQCD the shape function S has 
been calculated [55] ■ When all this is put together an ex- 
cellent description of data is achieved [63] , see Fig. [1] 



4 Production 

Production of heavy quarkonium is a far more complicated 
issue than spectroscopy and decay. This is due to the fact 
that in addition to the several scales which characterize 
the heavy quarkonium system, further scales due to the 
kinematics of the production process may also appear. 

4.1 Electromagnetic threshold production 

This is the simplest and best understood production pro- 
cess in the weak coupling regime, which is relevant for a 
precise measurement of top quark mass in the future In- 
ternational Linear Collider. The cross-section at NNLO is 
known for some time |65] and the log resummation at NLL 
is also available [33l[2T| . The effort is now in calculations 
at NNNLO, where partial results already exist [66 ^ 68 ( 167 ] 
[551[7ni[7T] . and in the log resummation at NNLL, where 
partial results also exist [68U72] . At this level of precision 
electroweak effect must also be taken into account [73 1 I74 ] 
[75] (see [751[77] for recent reviews). 

4.2 Inclusive production 

A factorization formula for the inclusive production of 
heavy quarkonium was put forward in [3], which was as- 
sumed to hold provided that the tranverse momentum pj^ 
was larger or of the order of the heavy quark mass, 

= ^i^(o|o^(M)|o), 

n Q 

The formula contains a partonic level cross section (F„ (^) / 
mg"""*) in which the heavy quark pair may be produced 
in color singlet or a color octet state, and long distance 
production matrix elements ((0|O^(/i)|0)) which encode 
the evolution of the pair to the actual physical state. The 
production matrix elements have the generic form 

o,f = x^/c„v ( ^ 5] Iff + + x| ) i^^icu 

\ X mj / 

where the sums are over the 2 J -f 1 spin states of the 
quarkonium H and over all other final-state particles X. 
K-n, K,'^ are gluonic operators (including no gluon con- 
tent and covariant derivatives). The production matrix 
elements were assigned sizes according to the velocity scal- 
ing rules of [3j, which, as discussed above, correspond to 
the weak coupling regime. A great success of this for- 
mula was the explanation of charmonium production a 
the Tevatron [78] , and it has been applied to a large num- 
ber of production processes (see [79l[80l[8T] for reviews), 
including some NLO calculations in photoproduction jSU 
[55] . This factorization formalism has received a closer look 
recently in the framework of fragmentation functions, and 
has been proved to be correct at NNLO in as{rnQ), pro- 
vided a slight redefinition of the matrix elements is car- 
ried out [M1[551[5B] . The interplay of the scales mqv, rngv'^ 
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and Aqcd has not been discussed in detail for production. 
This may be important in resolving the polarization puz- 
zle at the Tevatron [87], since the NRQCD prediction that 
heavy quarkonium must be produced transversely polar- 
ized in fragmentation processes |88j , sometimes phrased as 
a smoking gun for NRQCD, depends crucially on the ve- 
locity counting used. For instance, in the strong coupling 
regime one would not expect a sizable polarization 89j . 

Near certain kinematical end-points NRQCD produc- 
tion processes must be supplemented with coUinear de- 
grees of freedom, in analogy to semi-inclusive decays dis- 
cussed above. This issue has recently been addressed using 
SCET [QOllOTllMllMl. 



4.3 Exclusive production 

Although no detailed formalism has been worked out for 
exclusive production, the basic ideas of NRQCD factor- 
ization have indeed been applied to it f94','95] , mostly af- 
ter the surprisingly large double charmonium cross-section 
first measured at Belle (9^ . The traditional light-cone fac- 
torization does not seem to have major problems to ac- 
commodate this cross-section [97l!981l99]. although some 
modeling of the light-cone distribution amplitudes has 
so far been required. Recently, a further factorization of 
these objects in NRQCD has been presented jlOOj . and 
a detailed comparison of the NRQCD and light-cone ap- 
proaches has been carried out [lOlj (see also jl02j ). The 
latter suggests that the actual cross-section may well be 
accommodated into the NRQCD factorization results once 
higher order contributions in the velocity expansion are 
taken into account. 



5 Conclusions 

The NRQCD formalism provides a solid QCD-based frame- 
work where heavy quarkonium spectroscopy and inclu- 
sive decays can be systematically described starting from 
QCD. An NRQCD factorization formalism has also been 
put forward for semi-inclusive decays, inclusive produc- 
tion and, more recently, exclusive production, which, in 
spite of its successes, is not so well understood theoret- 
ically. Nevertheless, remarkable progress has been done 
recently concerning the structure of the factorization for- 
mulas for inclusive processes and the relation to the light- 
cone formalism for exclusive ones. Certain kinematical 
end-points both in semi-inclusive decays and production 
require the inclusion of coUinear degrees of freedom in 
NRQCD. Progress has also occurred here by combining 
NRQCD and SCET. 
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